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Abstract 

We regard the background of space-time as a physical system composed of discrete 
volume elements at the Planck scale and get the internal energy of space-time by Debye 
model. A temperature-dependent minimum energy limit of the particles is proposed from 
the thermal motion part of the internal energy. As decreases of the temperature caused 
by the expansion of the universe, more and more particles would be ”released” because of 
the change of the energy limit, we regard these new particles as a source of dark energy. 
The minimum energy limit also leads to a corrected number of particles in universe and 
a modified conservation equation. According to the modified conservation equation, an 


effective cosmological constant consistent with its observed value is obtained. 
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I. INTRODUCTION 


The astrophysical data obtained from high redshift surveys of supernovae [1, 2] 
reveals the fact that the expansion of the universe is accelerating. This acceleration 
is attributed to an unknown entity called dark energy (DE). In the past few decades, 
various modifications of general relativity have been proposed to explain the origin 
of it. Some of the most studied models are: the cosmological constant model, f(R) 
gravity [3-6], scalar-tensor theories [7, 8], quintessence models [9, 10] and phantom 
models [11, 12], etc. Unfortunately, as one of the most successful interpretations of 
DE, the value of the cosmological constant A calculated by the cosmological constant 
model is 60-120 orders of magnitude smaller than that given by quantum field the- 
ory [13]. This indicates that there is an irreconcilable contradiction between general 
relativity (GR) and quantum field theory (QFT). 

In order to eliminate the contradiction between these two theories, various quan- 
tum gravity theories [14-20] have been proposed in past few years. As an approach 
to find a suitable quantum gravitational theory, some physicists have considered the 
effect of the discreteness of space-time at a very small scale. In these theories, the 
fundamental description of space-time is only an approximation to manifold. For in- 
stance, causal set theory [21-24] described the space-time as a set of discrete points 
and developed the quantum dynamics of particles at these points. The points fixed 
on the background manifold are treated as the position of particles, therefore, the 
discreteness of the points will affect the trajectory of the particles at a tiny scale. A 
recent work [25] completed by Alejandro Perez and Daniel Sudarsky interpreted the 
cosmological as the accumulation of small violating the energy conservation which 
derived from the discreteness at the Planck scale. In these studies, the background 
of the universe was considered to have discrete microstructures, we follow this view 


and further hypothesize that the discrete points on the background have properties 


similar to particles, thus one can study the dynamics or thermodynamics of the back- 
ground. This assumption is equivalent to figure out the type of interaction between 
volume elements (voxels) of space-time, not just whether it is discrete or continuous. 
If there is a suitable micromodel to describe the background, it will provide possibil- 
ities for researching the energy exchange between background and the particles on 
the background, this is exactly the work studied in this article. 

In this paper, we regard the background as a thermodynamic system with strong 
interaction and discrete structure at the Planck scale, so the Debye model [26] can be 
used. By calculating the internal energy expression of the background under thermal 
equilibrium conditions, we obtain a temperature-dependent energy limit, considering 
the influence of the thermal motion of the background on the particles, this energy 
limit is assumed as the lowest energy of the particles. As the energy limit decreases, 
more and more particles contribute to gravity. As time goes on, this new energy will 
dominate the expansion of the universe. We regard it as a source of dark energy 
and calculate the cosmological constant corresponding to energy density of these 
new particles, which is consistent with the observation value of the cosmological 
constant [27]. 

The paper is organized as follows. Sec. II briefly review the method of obtaining 
the Debye model by canonical ensemble theory, then replace the ordinary atoms in 
solid with Planck-sized voxels and obtain the expression of the internal energy in the 
background of space-time. In Sec. III, the relationship between the number of the 
”real” particles and the total number of particles in universe is obtained using the 
lowest energy limit. Then we get a new conservation equation and estimated the 
value of the cosmological constant in Sec. IV. Finally Sec. V gives some conclusions 


and discussions. We will use natural units in which h = c = kp = 1 in this article. 


Il. THE DEBYE MODEL OF THE BACKGROUND OF THE UNIVERSE 


Physically, to take background of space-time as an object of research, we need 
an appropriate model to describe this system. Assuming space-time is discrete, 
meaning that space-time has internal structure, so we can study it with the theory 
of statistical physics. Intuitively, a system used to describe the background of the 
universe must not be a thermodynamic system with weak interactions like ideal gas, 
because such a system is not stable enough, the system will ”tear” when the energy 
is too high. Even a rest particle may also move chaotically while be embedded in 
such a background. Therefore, the background system should be dense, with strong 
interaction between discrete voxels. Fortunately, in statistical physics, there is an 
appropriate example to describe such a system, that is the Debye model. 

In general solids, there is strong interaction between the atoms, which makes 
the atoms have their own equilibrium positions and do micro-vibration near their 
equilibrium positions. Assuming that there are N atoms in the solid and each atom 
has 3 degrees of freedom, then the solid has 3N degrees of freedom. The total 
potential energy ¢ can be simplified as 
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where & is the displacement of the l-th degree of freedom from the equilibrium, and 
oo is the potential energy of the system when all atoms are in equilibrium. In addition 
to potential energy, each degree of freedom has corresponding microvibration energy. 


Therefore, the total energy of the system can be obtained 
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where Pe, /2m is the kinetic energy of l-th degree of freedom. In canonical coordinates, 


the above formula can be expressed as the sum of 3N independent simple harmonic 
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motions. As we all know, in quantum mechanics, the energy of such 3N independent 


simple harmonic motion can be expressed as 
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where, w; is the normal frequency and n; is quantum number describing the i-th 
simple harmonic motion. In this way, we get a form of energy that we can work on. 

In statistical physics, it is convenient to use the canonical ensemble theory to 
study systems with strong interactions. The partition function of the above system 


can be expressed as 
aaa 
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According to the canonical ensemble theory, the internal energy of the system can 


be obtained 
3N 
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Generally, ġo is negative, and its absolute value is greater than zero energy, so we get 
a negative binding energy Up which is independent of temperature. The second term 
on the right side of (5) represents the energy of the thermal motion of atoms. To get 
the result of internal energy specifically, we need to find the frequency spectrum of the 
simple harmonic motion. Debye regarded the solid as a continuous elastic medium 
and gave a spectrum of the system. The number of simple harmonic motions in the 


range of w to w + dw is 


D(w)dw = La (5 + 5) wdw, (7) 


where cı and cz represent the propagation velocity of longitudinal wave and trans- 
verse wave respectively. In addition, since such a system has only 3N simple harmonic 


vibrations, there must be an upper frequency limit, so the following formula holds 
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for this, we get 
uh = (9) 
where 
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As can be seen from (9), wp is related to the density of atoms and the velocity of the 
elastic wave, this fact is important to understand the relationship between scale and 
wp when we apply (9) to the background of univese. Using the Debye spectrum, (5) 


can be expressed as 
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In statistical physics, this integral has different results under different temperature. 


When the temperature of the system is much higher than the characteristic temper- 


ature 0p = wp, ie. T > Op, the result of (11) is 
U = Uo + 3NT. (12) 


At low temperatures T < p, the internal energy of the system is approximately 


equal to 
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the result in (13) indicates that the thermal motion part of internal energy is limited 


(13) 


by the characteristic temperature 0p at low temperature, this property guarantees 
the stability of strongly interacting systems. Similarly, stability can also be guaran- 


teed when we apply the above formula to the background of the universe. 


Now, we have reviewed the process to get the internal energy expression of 
the strongly interacting system in statistical physic. As mentioned above, we regard 
space-time as a physical object similar to a strongly interacting system, therefore,the 
expression in (11) will be used to represent the internal energy of the background. 
We will see that the thermal motion part of internal energy plays a crucial role in 


the production of dark energy. 


II. MINIMUM ENERGY LIMIT 


Let us study (11) from the view of cosmology. As mentioned above, the space- 
time is tread as a discrete system which has strong interaction between each small 
voxel like general solids, thus the internal energy expression in (11) can also be 
used to represent the internal energy of space-time. The first term on the right side 
of (11) Uo is the binding energy of the background of the universe. We assume that, 
in general, this part of the energy does not interact with the particles moving on the 
background, so it does not contribute to gravity. The last term of (11) represents the 
energy of the thermal motion of the voxels, T is the temperature when the background 
and particles are in thermal equilibrium. This means that there is energy exchange 
between the background and the particles. Assuming space-time is composed of N 


voxels, then the thermal motion energy of one voxel is 
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where @p is determined by the number density of voxels and the parameter B in (9). 
In this article, the scale of voxels is at the order of the Planck scale 4p, thus the 


corresponding number density n = T3, where T, is the Planck tempetature, then (9) 


T 


can be re-expressed as 
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we suppose that the speed of transverse wave cı and longitudinal waves c> is equal 
to the speed of light,i.e. c3 = co = c = 1. In fact, as long as cı and cy is in the 
order of c, the choice of the value of them has little effect on the result. As we can 
see from the above formula, 6p ~ Tp, this means that the relationship T < @p is 
always hold during most period of the universe, so (13) is the correct expression for 
the background of the universe. Substituting the result of (15) into (14), a more 
specific energy expression can be obtained 
n? Tt 
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which means that, on the cosmological scale, there is a vibration is in the order of 
T*/T3 in every point around the background. Different from the general case where 
only material exists, the vibration of voxels in background will have an important 
effect on the movement of particles because the particles are embedded in the space- 
time. Therefore, it is reasonable to believe that there is a minimum value of the 
energy of the particles given by the vibration of the background. The particles may 
not show up below this energy limit. 

In order to reflect the energy exchange caused by the minimum energy limit, it 
is useful to calculate the number of particles. The period is paid attention to during 
which the reheating of the electroweak transition is completed, the temperature of the 
universe is much higher than that of the static mass of most particles, where radiation 
is dominant. According to statistical physics, the state density of relativistic particles 
is given by 


D(e)de = ede, (17) 
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assuming that the particles follow the Boltzman distribution, in the presence of a 


minimum energy limit, the new number of particles is expressed as 
= YEc 
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where y is coupling constant that indicates the correlation of particles and space- 
time, and the value of y is around 1. €e is the characteristic energy obtained by (17). 
To get the result of the above formula, the expression of e7% must be known. It can 


be gotten from 
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Substituting the (17) into (19), one obtains 
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Then one gets the number of particles with energy above ye, from (18) is 

E _7T3 6 3 

Ñ=N-N |e Tý Ca aa (21) 
where Ñy = Ty, Planck temperature is the upper limit of temperature allowed, and 
even at the beginning of the time we studied, the temperature T, which approximates 
100 GeV, is much lower than the Planck temperature Tp ~ 10!° GeV. Ignoring the 


small terms in (21), a simplified modified number of particles is obtained 
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The above formula means that the ”real” particles in the universe are not the 
number of all particles (if there is a certain number of particles at the beginning of 
the Big Bang and it is roughly unchanged at a latter time), but partly cut off by a 
very small energy limit which changes with T’. As the temperature decreases, more 


and more particles are released, this process is accompanied by the accumulation of 


energy. When T — 0 and NWN , this energy accumulation process is over. In 
our study, this part of the energy is considered as the source of dark energy, and 
according to our calculations, we will see that the energy density of this part of 
energy is consistent with the energy density represented by cosmological constant 


observed. 


IV. DARK ENERGY 


At the end of the previous section, we obtained a corrected the number of 


particles. In this case, the corresponding particle number density is 
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In the standard cosmological model, the volume of universe is proportional to a(t)’, 
which means that the particle number density n œ a(t)”. Therefore, a corrected 
Hubble parameter is obtained 
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where H is the Hubble parameter in standard cosmological model. Now, we pay 
attention to the influence of such a Hubble parameter on the conservation of energy 


momentum tensor. Under the flat Friedmann-Lemaitre-Robertson-Walker metric 
ds? = —dt? + a(t)’ dz”, (25) 
the conservation equation V"T),, = 0 takes the form 
p+3d (1+w)p=0, (26) 


where, Tw is the energy momentum tensor of the perfect fluid, p refers to energy 
density (including radiation and matter) and w is the Equation-of-State (EoS) pa- 
rameter of the perfect fluid. Replacing the Hubble parameters H with (24), one gets 
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a corrected conservation equation as 
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Comparing with (26), a term that represents the energy flow density appears on the 
right side of the new conservation equation, defining 
= Paa 
a = -F7 (1 +w) zT P, (28) 
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on the large scale of the universe, the decrease of temperature leads to a negative T, 
if there is always —1 < w < i in every period of the universe, pa is always positive 
or close to zero. This means that this part of energy will gradually accumulate in 
some way and contribute to gravity, which we tend to treat it as the source of dark 
energy. 
Now, let us caculate the energy density accumulated by pa, then we can get the 
following integral 
to T2. To _ T2 
p=- f FUH) pioa =- f (1 +w) paged (29) 
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where, tren is the end time of the reheating period after the inflation caused by the 
electroweak transition, and tp is the time of the current universe. Ten and To refer 
to the thermal equilibrium temperature of the universe at the corresponding time 
respectively. At the time teh, the universe has just finished reheating, in which 
the rest mass of most particles is below the temperature T, en and the universe is 
dominated by radiation, the energy density is 

z2 


p= gT (30) 


where g, œ~ 100 is the effective degeneracy factor during the radiation-dominated 


period [28]. 
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At the time to, the current universe is dominated by the dark energy. Strictly, 
(29) is not applicable to the non-radiation-dominated universe, bacause (22) or (23) 
is obtained from the condition of relativistic particles. To get the correct number 
of particles in the non-radiation period, one just need to replace the state density 
in (18). For example, for non-relativistic particles, the state density is 


V 3/2 1/2 
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In this way, we will get a result slightly different from (29) and this result affected 
by temperature will be less than the period dominated by radiation. Today, since 
the Equation-of-State (EoS) parameter w — —1 and the energy density is much 
smaller than the energy density of the early universe, therefore, pa (To) in (29) can 
be ignored (although it is inaccurate for today’s universe). Thus the result of pa 
is actually determined by the temperature at the initial moment. Substituting (30) 


and the value of each constant into (29), one gets 
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therefore the cosmological constant A is 
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This result is consistent with Perez’s conclusion in [25], but from the completely 
different view. A slight difference to [25] is that the initial temperature we choose 
is the temperature after reheating, not the temperature Tew of the electroweak 
transition. As mentioned above, T,en is in the order of 100 GeV, but a more accurate 
value of T en is important for calculating the energy density of dark energy. Since 
T,eh is the temperature of the universe after reheating completed, its value should 
be slightly lower than the temperature Tw because of the expansion of the universe 


during reheating [29, 30]. According to our calculation, if the value of T,en is between 
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20 GeV and 30 GeV, we can get a result of A close to the observed value of the 


cosmological constant A,p5. 


V. CONCLUSIONS AND DISCUSSIONS 


In this work, we studied the background of the universe as a physical system 
and ulteriorly thought that the background has the microstructure at the Planck 
scale. Under the assumption that the background is composed of strongly interact- 
ing voxels, we used the Debye model and the canonical ensemble theory to describe 
the thermal motion of these voxels. With the consideration of the influence of back- 
ground on particle motion, we obtained a minimum energy limit of particles from the 
expression of internal energy. Just as the beach appears with the ebb tide of the sea, 
the particles with lower energy will gradually participate in the contribution to grav- 
ity as the energy limit drops. These newly emerged particles provide a new source 
of energy for the universe, which is referred to the dark energy. According to the 
relationship between the particle number density and the scale factor, we obtained a 
new conservation equation and calculated the energy density of the dark energy. The 
corresponding cosmological constant was also calculated which is consistent with the 
observation value of the cosmological constant. 

The discussion in this article shows that the energy density of dark energy is 
largely determined by the temperature at the end of the reheating period, which 
means that the dark energy may mainly comes from the accumulation of the early 
universe. As the decreases of the energy density of other components, dark energy 
would gradually dominate the expansion of the universe. Furthermore, such a result 
comes from a minimum energy limit caused by the microstructure of space-time, 


which reveals that space-time may have a discrete structure at the Planck scale. 
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